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We investigate the collision dynamics of two metastable strings which can be viewed as tube-like 
domain walls with winding numbers interpolating a false vacuum and a true vacuum. We find that 
depending on the relative angle and speed of two strings, instability of strings increases and the 
false vacuum is filled out by rapid expansion of the strings or of a remnant of the collision. 
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^ I. INTRODUCTION 

Qh It is becoming a growing consensus that the string landscape is one of the key concepts to reveal the early universe 

[U |2i . It suggests that vacuum structure of a field theory itself would be quite involved. Recently, inspired by this 

idea, a false vacuum has been widely exploited in phenomenological model building. In particular, there were drastic 

progresses on the supersymmetry breaking and its applications to phenomenology. (See [3j |4] for reviews.) Quite 

remarkably, in a wide class of models, a false vacuum is selected via a symmetry breaking associated with topological 

^^ soliton formation. A soliton is an energetic impurity in the false vacuum, which causes inhomogeneous false vacuum 

decay by a semiclassical process: Since the energy at the core of a soliton is high, if a lower vacuum exists, the 

^^ configuration of a field in the core may roll over the potential hill toward the lower vacuum and transform to tube-like 

D object. As the lower vacuum is energetically favored, the core of the tube-like soliton may start to expand immediately 

f^, 0^, depending on the parameters. In this case, the Universe is filled by the true vacuum via rapid expansion of the 

unstable soliton. 
^-H A cosmic string is one of the key solitons to probe information of the early universe or to prove (or exclude) a 

^ phenomenological model [10 . A fascinating aspect of the cosmic string is to make a string network. Short closed 
^^ strings are broken off from a long string and are constantly loosing the energy by, for example, gravitational wave 

or Nambu-Goldstone boson radiation. Then, it is believed that the energy fraction of the string finally reaches so- 
called the scaling regime. In making such a network, reconnection of strings that occurs frequently plays a crucial 
role. Therefore, it is important to investigate the dynamics of reconnection of (metastable) tube-like solitons. In this 
^ paper, we firstly show that a metastable string has the maximum winding number for stability, depending on the 

parameters of the potential. Then, we claim that even if solitons are metastable against rolling over the potential 
with an appropriate parameter choice, reconnection of such strings drastically increases instability and can induce 
inhomogeneous vacuum decay at the moment of the reconnection. We investigate stability/instability conditions for 
the strings under the reconnection by changing the relative velocity and the relative angle of two strings. 

Here, we focus on a simple single field model with a false vacuum and a true vacuum. This is partially because 
;J] the reconnection process is highly involved and is not easy to analyze complicated models with two or more fields by 
^ computer simulations in reliable precision. At any rate, it is very useful to reveal novel phenomena by exploiting a 
simple model which would be shared by a wide class of theories. 

The organization of this paper is as follows. In section II, we set up the model and show numerical solutions for 
metastable global strings. In section III, by using approximations, we estimate the maximum winding number of a 
metastable string and study analytically instability of colliding two strings. In section IV, we investigate the dynamics 
of colliding strings by simulations. We survey parameter dependence of instability by varying the angle and relative 
speed of the strings. Section V is devoted to conclusions and discussions. In Appendix A, we explain our scheme for 
numerical studies. 

II. SET-UP OF MODEL AND GLOBAL STRING 

To illustrate growing instability of metastable strings under reconnections and show generality of such phenomena, 
we consider a simple single field model with false and true vacua. This is an ideal example to demonstrate various 
aspects of the collision dynamics of two metastable solitons which would be common in a wide class of models. The 



Lagrangian of a complex scalar field X which carries a charge of global U(l) symmetry is given by 



c = \d^x\'' + v{x). (1) 



We engineer a false vacuum and a true vacuum by employing the following sixth order of the potential, 

V{X) = ^^'\X\^(5+^^{\X\^-r,'f 



(2) 



where J is a dimension- less constant determining the amplitude of a local minimum oiV{X) at X ^ 0. This potential 
has the global minimum at X = where F = 0, a local maximum at X = X^ax, and a local minimum at X = Xmin, 
which are given by 



K.i„ = V{X^^) = V*%^{C + 2)2(1 - c), K,ax = V(X„ax) = ^^^(c " 2f{c + 1), 
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l^minl = V\l^-, l^maxl = ^V ~3~' 



(3) 



where we introduced several dimension-less quantities, c = y^l — 3(5/C^, ( = r]/M and e = /i/r]. Figure 111 shows the 
schematic picture of this potential. The Euler-Lagrange equation we solve in section IV for simulations of colliding 
strings is given by 



^'^ AX + i^ = 0. 



df^ 



dX* 



(4) 




FIG. 1: The sixth order of potential Eq. ^. On the false vacuum at X — Xmin, the global U{1) symmetry is spontaneously 
broken. 



In the false vacuum, the global U{1) symmetry is spontaneously broken. We simply assume that a global string 
is formed. Consider a static cylindrical solution of the field equation in the cylindrical coordinate, (r, ^,z). First we 
decompose X into the radial and angular parts. 



X = r]R{r)e'''^ (n = 1, 2, 3, • • • )• 
Then the explicit form of the equation of R{r) is obtained from Eq. Q 



(5) 



(6) 



where x = errj. We look for the solution where the scalar field in the interior of the string stays in the true vacuum, 
and that stays in the meta-stable vacuum at the exterior. That is, such a solution should satisfy R = ai x = and 
R -^ Rmin at X -> oo. 

Here we emphasize that there can be (meta) stable tube-like solution. Since the string core has lower energy density 
than that of its exterior, larger string radius seems to be favorable for the system, which predicts the roll-over process. 
However, there must be the domain wall between its core and exterior, whose energy becomes larger for larger string 



radius. As a result, there arises a (meta)stable tube-like solution, depending on the parameters. We will see it in 
more detail in the next section. 

The static axial-symmetric solutions can be obtained by solving Eq. (IgI) with the successive over-relaxation method 
with the relaxation factor u = 1.0,0.5 and 0.3 for n = 1,2,3, respectively. See Sec. [A] for details of our numerical 
schemes. Using the boundary conditions i^(rb) = Rmin and R{0) = 0, we find the stable solutions. In Fig. [2J 
we plot the field configurations (left panel) and the potential energies (right panel) of the numerical solutions with 
e = 0.1, C = 4.0 and (5 = 1.0forn = l,2,3. We confirmed that the field configurations are insensitive to the position of 
boundary, rb, as long as it is sufficiently far from the domain wall. At the region where the potential energy becomes a 
peak, there is a domain wall which is the surface of the tube/cylinder. It is found that the radius of the tube/cylinder 
depends on the winding number, n, and the higher- winding solution tends to be thicker. Moreover, the thickness of 
the domain wall is quite insensitive to the winding number. 





FIG. 2: The radial field configuration R{r) (left) and their potential energies (right) with e = 0.1, (^ = 4.0, and 6 = 1.0. 



III. ANALYTIC ESTIMATIONS 



A. A schematic illustration of stability of the tube-like string and bubble 



In the previous section, we demonstrate numerical solutions for the static tube-like string. Before we consider the 
detailed estimation of their structure, we here give a rough approximation and clarify the stability of those solutions 
by using a simple thin-wall approximation. It is obvious that the stability of the tube-like string depends on the 
difference in the energy density of the true and false vacua 



Vf, 



"alse 



Vfr 



K, 



5ij:'7]\ 



(7) 



When e = ((5 = 0), the vacua inside and outside the string have the same energy. Then the tube-like string is stable 
because it is supported by the topological reason. On the other hand, when e is large enough, the core energy density 
wins out and the string grows and expands. Thus, one expects that there exists a metastable tube-like string in an 
intermediate range for e. Let us clarify this by using the thin-wall approximation where the wall of the string is thin 
compared to the core and the exterior. The thin wall's contribution to the energy density can be approximated by 
the constant surface energy, a. Then, the tension of the string with radius w is estimated by 



E(w) = n C log — + (jw 
w 



evo^ , 



(8) 



where the first term is the contribution of the fiux of the U{1) global current with Vc being a cutoff scale and C c^ 27rr]'^ 
being a constant of a mass dimension 2, the second term is of the wall and the last term is of the core of the tube-like 
string. Here we assume e, a, and C are independent of w. The dependence of the string energy on e is shown in Fig. [3] 
When e = 0, as we explained, there exists a local minimum, which implies that there exists a stable tube-like string. 



E{w) 





FIG. 3: A string energy (left) and a bubble energy (right) as functions of the core radius w. 



In a region < e < ec with Cc = 3§^, there exist a local minimum and a global maximum at 



Wn 



4g ' 



^max — 
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(9) 



In this case the tube-like strings exist but they are metastable. The tube-like string whose core radius is larger than 
'^max is unstable. When e exceeds the critical value, e > ec, there are no local minima at all. Namely, the string is 
unstable and the string core expands infinitely. In this paper the metastable strings are concerned. Although they are 
metastable as static configurations, they will become unstable in some dynamical processes such as string collision, 
annihilation and reconnections. 

From Eq. ([8| with n being zero, the fate of the 2 dimensional bubble can be also clarified. As shown in Fig. |3J the 
bubble never grows if e = 0. Once the positive e is turned on, the bubble becomes metastable. The critical radius 
of the bubble is Wc = ^. The bubble lager than Wc expands infinitely and the true vacua inside the core wins out. 
There are no metastable tube-like solutions. 

The thin-wall approximation explained above is simple and gives an insight about the stability problem. Nev- 
ertheless, it is of only limited accuracy and we need a better analysis to get more quantitative results beyond the 
qualitative properties. To go beyond the thin- wall approximation, we will work out another analytical study in the 
next subsections. 



B. thickness of walls 



Now let us go beyond the thin- wall approximation and consider the stability of the tube-like solution in more detail. 
From the numerical calculations of stable static solutions, we found that the tube-like solution is well characterized 
by the radius of cylinders, w^ the winding number, and the thickness of the walls, d. Since the thickness of the wall 
is quite insensitive to other two parameters, it is possible to estimate d from the variation of the total energy E with 
respect to (i, and it can be done by considering only the case of n = 1. 

We here simply approximate the solution with n = 1 as the following piecewise linear function 




r < d 
d <r 



(10) 



where we set i^ = 0, see Fig. [2] Then, we approximate the volume integral of the potential energy and gradient 
energies as 

-^potential — ^^ V »^max »^minJ5 
-LVprradr TT/y -''^niin ' 



^grad,r 



^grad,e = 7rr?2i?^i„(l + 21og^) 



(11) 

(12) 
(13) 



The derivative of the summation of these energy components with respect to d gives the desired value of d. 
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(14) 



C. upper bound of winding number 



We consider how the radius of the cyhnder is determined and the upper Umit of the winding number with which 
the static solution exists. As shown in Fig.|2J for n > 2, there are two parameters characterizing a cyhnder, the radius 
w and the thickness d of the surface. However, as shown in the previous subsections, the latter one is not sensitive to 
the winding number and is determined from Fmax — "Knin- Hence we here focus on the radius w of the cylinder while 
the width of the wall d is assumed to be the same as one given in Eq. (14). 

We simply approximate the solution and the potential as 



R{r) ^ < 





r — w 

d 



Rr, 



.V{r)^{ 







r < w 

^max W <r <W 

"141111 w^d<r 



(15) 



The volume integral of the potential energy is approximated as 



where we considered only the deviation from Vm 
into the radial and angular parts. 



^potential = -"KninTT^^ + 7rd{2w + C^)(14iax " "Knin), (16) 

On the other hand, that of the gradient energy can be divided 



^grad,r = 2^7/' / R\r)^ vdv = TTT]^ R^,^ [2^ + 1 
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where Vc is the cut-off length, Tc ^ oo, and it should be properly regularized. Then the total energy becomes 



E = E, 



potential 



E, 



grad,r 



■E, 



grad,^- 



(19) 



Let us look for the values of w to minimize the total energy. One can easily see that there is no global minimum of 



E since the term being proportional to —w in Eq. (16) implies that E -^ — oo for w ^ oo. Instead, we investigate 



whether there is a local minimum in the region, w > ^. The derivative of E with respect to w is calculated as 

2^rfRl 



dE 

dw 



•"mm 



G{w), 



where 



G{w) = -2p- +2 + 2n 



d'Vrr. 



-1 



log 1 



(20) 
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What we have to do is to find w to satisfy G(i^) = for i^ > 0. We plot E{w) and G{w) with specific values of e, (" and 
6 in Fig. [4J From the right panel of this figure, it is found that there are two zero-points for n < 6, and the smaller 
one is the desired value of k;, at which E is locally minimized, and the other zero-point gives the local maximum of 
E (see the corresponding lines in the left panel). Note that the local minimum of n = 1 is k; = 0, which is consistent 
with the numerical solution in Fig. [2] 

In order for these zero-points to exist, the local maximum of G{w) should be positive. Actually, G{w) for n = 7 is 
always negative, and in the left panel of Fig.plthe line for E{w) with n = 7 has no local minimum, which indicates the 
cylinder solution with n = 7 is unstable. To proceed this calculation, we assume d/w <C 1 to expand the logarithm 
function. Then we find the G{w) has a local maximum. 



G{wc) - 2 - 2nv/2^, 
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Therefore the condition on the winding number required for the local minimum of E is 
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FIG. 4: An example of E(w) / (2iiTf R^^^) {left) and G{w) [right) with specific values of e = 0.1, C = 4.0, and 5 — 1.0. 



D. Instability of spherical bubble and critical volume 



In subsection III A we have seen that the two-dimensional bubble is unstable and it shrinks (expands) when its radius 
is smaller (larger) than the critical value. Here, we study the spherically-symmetric configuration (three-dimensional 
bubble) of the scalar field separated by the domain wall. It is also expected to be always unstable. After two cylinders 
collide with each other, there appears a spherical object at the impact point. To study the instability of such an 
object, we consider the following ansatz in the spherical coordinate, 
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with assuming that d is given by Eq. (14). We simply assume that the field is homogeneous in the azimuthal and the 



polar directions^. The potential and radial gradient energy are approximated as 



E. 



potential 



rTTW^V^in + :^Tt{{w + df - W^KKiax " Klin), 
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{{w + d)-^ -W^}ri^, 



(25) 
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and thus the total energy becomes 



E{W) = -TT [-W^V^in + 2{{W + df - W^}{V^^^ - Knin)] , 



(27) 



where we used Eq. (14) to eliminate i^min- The function E{w) has the only peak dit w = Wc in w > 0, and the critical 
radius Wc is given by solving dE/dw = 0, 



Wc 



1 + vTTp 



d, 



(28) 



where p is defined in Eq. (21). For w < Wc the radius of the sphere starts to shrink since dE/dw > 0. On the other 
hand, for w > Wc the radius goes to the positive infinity since dE/dw < 0. That is, if the volume of the spherical 
object is larger than the critical volume Airw^/S, this grows infinitely. 



^ Strictly speaking, the separation of variables is not justified, since it is impossible to expand X in the spherical harmonics due to the 
non-linear terms in the potential. 
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FIG. 5: The volume colliding with another cylinder. 



Now we give an analytical estimate of the (in) stability of the colliding cylinders. Let us consider the situation that 
two cylinders collide with the collision angle a as shown in Fig. [5] Assuming that a spherical object is produced via 
the fusion of two cylinders after collision and its volume is equal to the total volume of the two colliding segments of 
cylinders, 2l^eg7 the initial volume of the spherical object is calculated as 



^otal = 2T4eg = 

Then the condition for the sphere to grow infinitely is given by 

47r 



Attw^ 
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-wl < Vt 



total 



sin a < 3 
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where w is obtained by solving G{w) = in Eq. (20), and hence assuming d/w <C 1, this is calculated as 



1- a/1 - 2n^p , 

w = a. 

2p 



(29) 
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(31) 



As a result, we obtain the upper limit of the collision angle so that the spherical object created at the impact point 
can grow infinitely. 



3/1- a/1 - 2n2p 

sma < - — , — 

8 I 1 + vTTp 



(32) 



For example, e = 0.1, C = 4.0 and 5 = 1.0 give the upper limit, a < 7.78 x 10 ^, for n = 2, and a < 2.50 x 10 ^ for 
n = 6. 



IV. SIMULATION SETUP 



We explore the (in)stability of collision processes of two meta-stable strings. The individual string is given by 
the axial-symmetric solution which is numerically obtained by solving Eq. Q without the time-derivative term, 
corresponding to Eq. (|6|. This procedure is nothing but the one-dimensional boundary value problem, and is carried 
out with the Gauss-Seidel method. Each of the resultant strings is put on the x-z surface parallel to each other 
separated a bit. The strings are Lorentz-boosted to collide with each other at an angle with a measured on the 



x-z surface including the origin (to be the impact point when colhding). The free model parameters are e and (", 
which control the strength of the self-coupling and the energy scale of the phase transition, respectively. The fiducial 
parameters are listed in Table [T[ Note that we choose ( = 2.2 for numerical simulations, while (" = 4.0 has been used 
so far, to enhance the instability. 

The initial separation between the strings is fixed to be 20?7~^, while the radius of the string is ~ 5r]~^ for the 
fiducial choice of parameters. According to Ref. [10], the superposition of the two strings is given by 



X = 



XiXo 



r]Rn 



(33) 



where the denominator is determined by the dimension analysis and the fact that |X| ^ TjRmin far away from two 
strings. 

Then we solve Eq. Q in the three-dimensional cartesian coordinate with the Neumann conditions on the boundaries, 
9X/9n* I boundary = 0, whcrc 71* is the normal vector to each boundary. We use the Leap-frog method for the time 
domain and approximate the spatial derivatives by the 2nd-order central finite difference. The simulations are stopped 
at t = L/2, when the information at the impact point arrives at the nearest boundary of the computational domain. 



model parameters 


simulation parameters 


e 


0.1 


grid size 


N 


480^ 


c 


2.2 


box size 


L 


240ry-^ 


6 


1.0 


time interval 


At 


0.2ry-^ 


ni 


1 


total steps 


Ns 


600 


712 


1 


simulation time 


tf 


120ry-^ 



TABLE L The fiducial parameters. 



V. RESULTS 



Through the numerical simulations, we find that the collision processes end up with either simply reconnecting 
and going away from each other, or creating unstable objects with a higher winding number at the impact point. 



Figs. [6pT] show the snapshots of the computational domain during simulations, the leftmost panel is at the initial 
time t = 0, and the rightmost one at the final time t = tf. The surfaces in them represent X = X^ax given in 
Eq. (|3|, and hence in the interior of them the field X lies in the true vacuum. The first four figures. Figs. [6][9J are 
the cases of parallel strings with relatively small a, and the remaining ones. Figs. 10 and 11 , are those of anti-parallel 
strings with a ^ tt. Moreover, we plot the field configurations and those phases on the surface at z = L/2 and 
t = tf in Figs. 12 and 13 In the unstable cases, PI, P3, P4 and Al, we find the true vacuum homogeneously spreads 
over the interior of thebubbles. As for the phases, one can observe that the winding number is conserved. For the 
parallel pairs, one can find that the total winding number isn = l + l = 2 (two sets of blue^green^red^blue), 
if one follows the trajectory around the bubble. In the interior, there seems to be a large number of points around 
which the phase is rotated, although the winding vanishes if one follows the trajectory around the impact point. On 
the other hand, for the anti-parallel pairs, the total winding number around the bubble becomes n = 1 — 1 = 
(green^blue^green^red^blue^red^green.) 




FIG. 6: [PI] Failed reconnection of parallel string pair with {v,a) = (0.82c, 0.027r). See [11 for full animations. 



The (in)stability of a resultant string after collision depends on the velocity of two cylinders, v^ and the collision 
angle between them, a. We surveyed the parameter space (a, v) to check the stability. Fig. 14 shows the stability of 




FIG. 7: [P2] Successful reconnection of parallel string pair with {v,a) = (0.90c, 0.027r). See [TT] for full animations. 




FIG. 8: [P3] Failed reconnection of parallel string pair with {v,a) = (0.94c, 0.027r). See [11 for full animations. 




FIG. 9: [P4] Failed reconnection of parallel string pair with {v,a) = (0.95c, O.IOtt). See [11 for full animations. 




FIG. 10: [Al] Anti-parallel string pair with {v,a) = (0.62c, 0.947r). Most of segments are annihilated. See JJj for full 
animations. 



parallel string pairs, and FigjTSJthat of anti-parallel pairs where a is close to tt. In order to systematically judge the 
resultant stability, we calculate two quantities and set a criterion for each. One is the growth rate of the volume of 
the true vacuum region during simulation. 
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and the other is the ratio of time when the true vacuum region grows to the total simulation time. 
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FIG. 11: [A2] Anti-parallel string pair with {v,a) = (0.90c, 0.947r). After collision, the true vacuum region is created around 
the impact point, and it starts to grow exponentially. See pj] for full animations. 




FIG. 12: The field configurations of [P1]-[P4] on z = L/2 (upper panels), and those of [Al], [A2] on y = L/2(lower panels) at 
t — tf. The colour contour represents \X\/ Rmin- The white colour indicates the overshoot, |-^| > i^min- 




0.25 0.5 0.75 



FIG. 13: The phases of [P1]-[P4] on 2; = L/2 (upper panels), and those of [Al], [A2] on 2/ = L/2 (lower panels) at t = t/. 



where ft returns the number of elements, At is the time interval, and tm represents the discrete time, tm = mAt. If 
and only if both /^i > tzic and hZ2 > ^2c are satisfied, we decide that the resultant strings become unstable, and thus 



11 

the true vacuum region grows exponentially. We set kic = 10 and hZ2c = 0-8- The second quantity, hZ2^ measures how 
monotonic the instability grows. After all, however, most of (in) stabilities can be determined by the first one, Hii. 

In Figs. [M] and [ISJ we plot red filled circles for unstable pairs, and green crosses for stable pairs. We find that the 
slow collision with small a results in the failure of the reconnection and thus being unstable, and quite high-speed 
collision can also make the system unstable in both cases with small (parallel pair) and large a (anti-parallel pair). In 
these figures, we also plot black circles labelled as PI, 2, 3, 4 and Al,2. The corresponding figures were already shown 
in Figs. [6pT| which are the representatives of each stable or unstable region in Figs. [M] and [151 

We explain in detail the mechanism of stable/unstable reconnection processes using Figs. 6pT Firstly, for the 
parallel pairs, our findings are followings. 



PI (relatively slow collision) : Two strings merge at the instance of impact and cannot be separated from each 
other due to the small collision velocity. Then the effective volume of the merged strings around the impact 
point exceeds the critical volume mentioned in Sec|III[ thus making the system unstable. 



P2 (moderate-speed collision) : Two strings can safely reconnect with each other. As shown in the middle panel and 
its right neighbour, a small bubble is created at the impact points of the strings. However, since the reconnected 
strings rapidly go away from the bubble and the volume of the bubble is relatively small, the bubble finally 
shrinks. As a result, the system becomes stable. 

P3 (high-speed collision) : Due to the large collision velocity, the kinetic energy induced to the bubble at the impact 
point is large enough to infiate the bubble until its volume exceeds the critical one. Finally the system becomes 
unstable. 

P4 (large angle high-speed collision ) : Since the collision angle is large, the initial impact of strings cannot create 
a larger bubble at the impact point than the critical volume. However, induced kinetic energy to the bubble is 
too large, and then the bubble can rapidly expand (the third panel). At the same time, the reconnected strings 
also rapidly go away from the impact point. As a result, the bubble is about to break up (forth panel). Finally, 
the expansion rate of the bubble overcomes (last panel). 

As for the anti-parallel pairs, they annihilate with each other in most of cases due to the cancellation of their 
winding numbers [Al]. It is, however, surprising that the high-speed collision makes the system unstable. Actually, 
after the collision, a bubble with no windings is excited at the impact point, and it grows exponentially. 
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FIG. 14: The parameter region for stable/unstable collision process for string-string collision. 



VI. CONCLUSION 



We have studied the dynamical instabilities of metastable strings after their collision in a model where the potential 
of a complex scalar field has a false vacuum state. 
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FIG. 15: The parameter region for stable/unstable collision process for string- antistring collision. 



Before performing numerical simulations, we analytically investigated in a simplified model the thickness of the 
domain wall constituting the surface of strings, and the existence of static solutions with large winding numbers from 
the viewpoint of energetics. As a result, we found that the thickness is determined from the shape of given potential 
and independent to the winding number, and there exists the upper bound of winding number for the stable solutions. 
Furthermore, in the same manner, we investigated the collision process of two metastable strings. Then we found that 
there is the critical volume so that the instability grows, and thus in some cases with a small collision angle and/or 
small velocities of strings, the reconnection is failed and a true vacuum bubble created at the impact point starts to 
grow infinitely. 

After that, we performed the three-dimensional field-theoretic simulation of two colliding metastable strings in 
the cartesian. The initial condition is given as the superposition of two Lorentz-boosted strings obtained in another 
numerical way. We surveyed the (in) stabilities of the collision processes on the parameter space (v,a), where v is the 
velocity of strings and a is the collision angle. 

Consequently, we found that the instability cannot always be observed for the string-string pairs. We fixed the 
winding number of strings as n = 1 and tuned the parameters controlling the shape of potential, e, (" and 5^ so 
that the static solution with n = 2 does not exist. Nevertheless the strings with most of combinations (v, a) are 
successfully reconnected, or, in other words, the parameter region in the (v, a) space where the instability grows is 
highly restricted. This fact indicates that the naive expectation, such that the instability would grow due to the 
temporal formation of unstable n = 2 strings at the moment of impact, is not true, and instead, our analytical studies 
before simulations, showing the existence of the critical volume to grow instability and the winding number itself is 
not responsible for the (in) stability, is well supported. In reality, although the unstable region non-trivially spreads 
over the (v, a) space as shown in Fig. 14, all of them can be explained by the critical volume. 

For the string-antistring pairs, the collision processes were basically stable, and the string annihilated with the 
antistring as expected. However, even in such system we found the unstable region for large v as shown in Fig. 15 It 



is surprising, but would be explained by the shorter time scale of the expansion of the zero- winding bubble created 
at the impact point, which eventually be larger than the critical volume, than that of annihilation. 
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Appendix A: numerical schemes 

To find the static configuration of an axially- symmetric metatable string, we solve Eq. (|6| with the succesive 
over-relaxation method. 

First of all, we discretise the coordinate such as Xj = jAx for j = 0, 1, . . . , A/", and represent the discretised R{x) 
as a vector Rj = R{xj). The spatial derivatives in the left-hand side of Eq. (|6| is replaced by the corresponding 
2nd-order finite differences. 



d^R 
dx'^ 



Rj-\-i ~ 2iij + Rj—1 dR 

Ax^ ' dx 



^3^^ -^3-1 (Al) 



We denote the other terms depending on R{x) in Eq. (|6| by S[R{x)] and those evaluated at x = Xj by Sj [R\ = S[R{xj) 
Then the equation to be solved becomes 



Rj-\-\ — 2^Rj H~ Rj—i /tj_|_i — Rj—i 

Ax^ 2xjAx 



= Sj[R]. (A2) 



Solving this equations with respect to Rj in the first term, we find trial values of Rj denoted by i?* , 

Rj±l^rRj-l A^^ f Q rol ^j + i - Rj-l 



Notice that Rj in Sj[R] is not R*. On each stage of iterations, we update Rj so that 

i?5"+'^ =ojR*j + {l- uj)Rf^ . (A4) 

This scheme with cj = 1 is equivalent to the Gauss-Seidel method. When we solve Eq. (p| for (^ = 4.0 in Sec. Bl| we 
set cj = 1, 0.5 and 0.3 for n = 1, 2, 3, respectively. For C, = 2.2 in Sec. [TVl we set a; = 1. 
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